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SUMMARY 


In an effort to improve understanding of the effects of variations in 
background water on reflectance spectra, laboratory tests have been conducted 
with various concentrations of sewage sludge diluted with several types of 
background water. The results from these tests indicate that reflectance 
spectra for sewage-sludge mixtures are dependent upon the reflectance of the 
background water. Both the ratio of sewage-sludge reflectance to background- 
water reflectance and the ratio of the difference in reflectance to background- 
water reflectance show spectral variation for different turbid background 
waters. The difference in reflectance is the only parameter considered in 
this study which does not have major spectral variations with different 
background waters. 


INTRODUCTION 

The National Aeronautics and Space Administration (NASA) , in cooperation 
with the Environmental Protection Agency (EPA) and the National Oceanic and 
Atmospheric Administration (NOAA) , is conducting a research program to evaluate 
the feasibility of remotely monitoring ocean dumping of waste products such as 
industrial wastes and sewage sludge. The program consists of both field exper- 
iments (refs. 1 to 8) and laboratory studies (refs. 8 to 10). As a result of 
the expense involved, the field results have generally been limited to one 
experiment with a specific discharge. Repeat experiments to investigate the 
reliability of data-analysis algorithms and the effects of variations in back- 
ground water and atmospheric interference have not been conducted. In an 
effort to improve this situation, a series of laboratory tests has been con- 
ducted to investigate the effects of variations in background water on a single 
pollutant, dilute pr imary-treated sewage sludge from the Southwest Philadelphia 
Waste Treatment Facility. It is the purpose of this paper to describe how 
upwelled reflectance spectra of dilute sewage-sludge mixtures changed as the 
content of chlorophyll a and suspended solids in the background water varied. 

In addition, several types of data-analysis algorithms are applied to the labo- 
ratory data to investigate the sensitivity of each algorithm to variations in 
background water. It is expected that such information from laboratory tests 
will aid in the selection of algorithms to be used in future field experiments. 

Use of trade names or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 


A 


SYMBOLS 

area of spectrometer entrance slit, cm^ 
[pu ( 1 bw^so ^ Ibw' iKW/cm^-sr-nm 


B 


c 


La(X)/Ta{X) , rtW/cm^-sr-nm 

vertical displacement of oscilloscope measurement, cm 


D 

E(X) spectral ir radiance, W/m^-nm 

K ratio of instrument throughput to vertical-scale sensitivity factor, 

cm^-sr 
mW/nm-cm 

La(X) upwelled radiance from light scattered into receiver by the 

atmosphere, nW/cm^-sr-nm 

Lrd(X) upwelled radiance from specular reflection of diffuse skylight, 
mW/cm^-sr-nm 

Li-s(X) upwelled radiance frcxn specular reflection of sunlight, mW/cm^-sr-nm 

Lgo(X) upwelled radiance above water surface from 100-percent diffuse 

reflector, mW/cm^-sr-nm 

Lu(X) inherent upwelled spectral radiance above water surface, mW/cm^-sr-nm 

L 2 (X) apparent upwelled spectral radiance at altitude z, mW/cm^-sr-nm 

A[L 2 (X)]gg = [L 2 (X)]gg - [L 2 (X)]j 3 ^,, iiiW/cm^-sr-nm 
P(X) spectral power, mW/nm 

S vertical- scale sensitivity factor, mW/cm-nm 

Ta(X) atmospheric transmission fraction 

a beam attenuation coefficient at X = 51 0 run, m“^ 

X wavelength, run 

Py(X) inherent upwelled spectral reflectance above water surface, percent 

of input (relative to a 100-percent diffuse reflector) 

A[Pu(X)lss = tPu(^)^ss - fPu(^)^bwf percent of input 

A[pu(X)]gg = [Pu(X)]gg ~ [ P u ( ^ ) 1 ss =4ppm ' percent of input 

a standard deviation of measurement error 

acceptance solid angle of spectrometer, sr 


2 



Subscripts: 


bw background water 

ss sewage sludge mixed with background water 


EXPERIMENTAL METHOD 

The arrangement shown in figure 1 was used for the laboratory tests. 

Major parts of the system include a water tank, circulation system, filtration 
and deionization system, solar simulator, first-surface mirror, and rapid-scan 
spectrometer. The light source is a 2.5-kW xenon short-arc lamp which produces 
a spectrum similar to that of the Sun at a solar elevation angle of 30°, as 
shown in figure 2. A more complete description of the laboratory and equipment 
is given in appendix A. 

Approximately 4 weeks prior to the laboratory spectral-signature tests, 
cultures of Phaeodactylon tr icornutum (a golden-brown algae) were added to a 
2-liter laboratory container of artificial seawater which had been prepared in 
accordance with reference 11. (See table I.) The nutrients in table II were 
added and artificial sunlight was used to grow the algae in a laboratory envi- 
ronment. When cell count increased such that chlorophyll a concentration 
exceeded 500 Vig/il, the mixture was transferred to a 95-liter tank which had 
previously been filled with artificial seawater. Nutrients were added to the 
diluted algae mixture and growth was continued under artificial sunlight for 
several weeks. When chlorophyll a concentration was again in excess of 
500 Ug/Si, the mixture was transferred to a 341 -liter tank and growth was 
allowed to continue. The desire was to achieve chlorophyll a concentrations 
between 700 and 800 yg/Jl in the 341-liter container. Portions of this mixture 
would then be added to the 11 600-liter tank shown in figure 1 to achieve mix- 
tures with the desired algae content. 

The 11 600-liter tank was then filled with tap water that had been fil- 
tered and deionized such that it contained less than 0.5 ppm of total suspended 
solids and less than 2 ppm of dissolved ionic substances. The chemical com- 
pounds in table I were added in sequence in the amounts shown after each pre- 
ceding compound had been observed to dissolve in the water mixture. The final 
artificial seawater mixture was then filtered for 15 hours such that concentra- 
tion of total suspended solids was less than 0.1 ppm and concentration of chlo- 
rophyll a was less than 0.2 Vg/H. 

The sequence of laboratory tests is summarized in figure 3. The high 
value of total suspended solids shown for artificial seawater plus nutrients 
(mixture 0) is believed to have been caused by incomplete dissolution of the 
powdered nutrients in the artificial seawater prior to the addition of algae. 

A portion of the mixture in the 341 -liter tank (chlorophyll a = 740 yg/J-, and 
total suspended solids = 43 pizn) was added to the 11 600- liter tank. The 
resulting water had a chlorojAiyll a concentration of 7 yg/S- and is labeled mix- 
ture 1 in figure 3 and table III. This is the first mixture for which upwelled 



spectral-reflectance measurements are presented. Two additional water trans- 
fers were made later in the day to obtain mixtures 2 and 3, which had chloro- 
phyll a concentrations of 15 and 20 yg/il, respectively. Artificial light was 
then added and the algae allowed to grow for 17 hours to obtain mixture 4. The 
next morning, analyses were made of mixture 4, which had a chlorophyll a con- 
centration of 31 Vg/H and a total suspended- solids concentration of 11.6 ppm. 
After measurement of the upwelled reflectance spectra, two additions of primary- 
treated sewage sludge were made such that total suspended-solids concentration 
increased frcrni 11.6 ppm in mixture 4 to 19.3 ppm in mixture 5 and 48.7 ppm in 
mixture 6. Table III indicates that algae cell count remained essentially 
constant . 

Upwelled radiance over the wavelength range from 400 to 980 nm was mea- 
sured for each of the six water mixtures in table III over the 2-day period 
of the tests. For each water mixture, diffuse reflectance measurements of the 
solar-simulator spectral input at the water surface were made with a 99-percent- 
reflectance white card to monitor optical stability of the laboratory system 
and to derive spectral reflectance from upwelled radiance values as described 
in appendix A. 


RESULTS AND DISCUSSION 
Artificial Seawater Tests 

Measured spectral radiance and reflectance curves are shown in fig- 
ures 4(a) and 4(b), respectively, for mixtures 1 through 4. All data are for a 
spectral resolution equal to 32 nm. Estimated values for standard deviation O 
are based on manufacturer's specifications (see appendix A) combined with data 
readup inaccuracy. Not included in the estimates of standard deviation are the 
uncertainties in the concentrations of chlorophyll £ and total suspended solids 
inherent in laboratory analyses of water samples. Standard deviation of chlo- 
rophyll a concentrations is estimated at 20 percent of the observed value. 

(All chlorophyll a values given in this report are the average of three mea- 
surements.) Estimated standard deviation of total suspended-solids values is 
25 percent at 7 ppm and 10 percent at 50 ppm. Variations in the spectral- 
radiance curves (fig. 4(a)) are caused in part by the spectral characteristics 
of the xenon-lamp input source (see fig. 2(b)). Figure 4(b) shows spectral- 
reflectance values which cancel the effects of the lamp spectra. Figure 4(b) 
shows little change in reflectance for chlorophyll £ concentrations between 7 
and 20 Ug/S- for the Phaeodactylon tricornutum algae species. As higher con- 
centrations of chlorophyll a are approached with mixture 4, a larger increase 
in spectral reflectance is observed. 

Spectral- radiance and reflectance values are shown in figures 5(a) 
and 5(b), respectively, for mixtures 4, 5, and 6. Mixture 4 is considered to 
be a turbid background-water mixture to which two quantities of sewage sludge 
were added to obtain mixtures 5 and 6. In figure 5(b), it appears that the 
addition of sewage sludge from the Southwest Philadelphia Waste Treatment 
Facility causes reflectance to increase by a nearly constant amount, within the 
accuracy of the measurements over the spectral range from 400 to 980 nm. 
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Ccjmpar ison With Previous Results 

Dilute pr im ary- treated sewage sludge from the Southwest Philadelphia Waste 
Treatment Facility has previously been tested and the results reported in ref- 
erence 9. Those tests were conducted with both filtered, deionized tap water 
and natural water pumped frcan the Back River at Hampton, Virginia. Figure 6(a) 
shows reflectance spectra from the reference 9 tests compared with the curve 
from figure 5(b) for the mixture with a sewage-sludge concentration of 37 ppm 
(mixture 6) . This curve for mixture 6 will be considered as representative of 
a 34-ppm mixture for the remainder of this publication, since the difference is 
within repeatability of water sample analysis. Figure 6(a) indicates that 
reflectance spectra for the two samples of sewage sludge at a concentration of 
34 ppm are different. Figure 6(b) shows reflectance spectra for two of the 
background waters, mixture 4 and Back River water. A reflectance spectrum for 
filtered, deionized tap water is not shown because of limitations of the labo- 
ratory system. (See appendix A.) Ccanparison of figure 6(a) with figure 6(b) 
indicates that the spectral-reflectance variation of each sewage-sludge mixture 
is quite similar to that of the background water for both mixture 4 and Back 
River water. 


Application of Selected Algorithms 

As discussed previously, one purpose of this study was to investigate the 
sensitivity to variations in background water of several algorithms for analy- 
sis of field data by use of laboratory reflectance data. In order to accom- 
plish this, it is first necessary to relate t 5 >welled radiance measurements at 
altitude to surface-reflectance values. An analytical expression for this 
relation is given in appendix B, along with similar equations for several 
algorithms. In the following paragraphs, laboratory reflectance values will 
be applied to these algorithms in an effort to infer probable characteristics 
of remote-sensing radiance measurements from aircraft or satellites. 

The ratio of sewage-sludge to background-water apparent radiance is 
related to the ratio of sewage-sludge reflectance to background- water inherent 
reflectance as shown in the following equation (which is equal to eq. (B3) ) : 

fl*z ( ^ ss ^ ^ Pu ( ^ ss/ ^Pu ( ^ bw^® ^ 

1 T“i — = 

iLz(X) ]5„ B + C 

Application of laboratory values from figure 6 produces the spectra of inherent- 
reflectance ratios shown in figure 7. Figure 7(a) shows curves based on contin- 
uous spectrometer data. Simulated values for bands 2 through 9 of the Modular 
Multispectral Scanner (M^s, Bendix Corp.) instrument have also been calculated 
and are shown in figure 7(b). These calculations were made by integrating the 
values of figure 7(a) over the bandwidths of 40 to 100 run for each channel of 
the M^S instrument. The M^S instrument was chosen for simulation because of 
its wide use in ocean-dumping remote-sensing experiments (refs. 5 to 7, for 
example). Figure 7(a) indicates a difference in [Pu (X) 1 gs/ fPu ^ bw 
different turbid background waters. The sharp peaks observed in the laboratory 
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spectra are not evident in the simulated M^s spectra because of the poorer 
spectral resolution of that instrument. If the spectral characteristics of 
the ratio of sewage-sludge reflectance to background -water reflectance vary 
with differences in turbid background water, it is reasonable to expect that 
the ratio of sewage-sludge radiance to background-water apparent radiance at 
altitude z will also vary according to equation (1). 

Another algorithm which has some potential for identification of marine 
pollutants is one in which the difference between the sewage-sludge radiance 
and background-water apparent radiance A[L 2 (X)]gs is expressed as a ratio 
to the apparent radiance of the background water [L 2 (X)]| 3 „. This ratio is 
related to inherent-reflectance values by the following expression (which is 
the same as eg. (B6)); 

A [ (X)]gg {A[P(j(X)l gs/ f Pu (X ) 1 

7 — n — = ^2) 

1 (X ) ] ^ ^ C 

B 

The reflectance ratio A[py(X) ] gg/[py(X) ] shown in figure 8 has been com- 
puted frcsn values in figure 6. From these curves, it is concluded that differ- 
ences in both reflectance magnitude and spectral variation occur for different 
turbid background waters. It is reasonable to expect that the ratio of the 
difference in apparent radiance to the background-water apparent radiance will 
also vary for turbid background waters. 

The third algorithm investigated was the parameter for the difference in 
apparent radiance A[L 2 (X)]gg. This parameter is related to the difference 
in inherent reflectance by the following expression (which is the same as 
eg. (B9) ) : 


A[L2{X)]gs = [Ta(X)] [Lgo(X)]{A[pu(X)]gg} (3) 

Again by use of values from figure 6, the inherent-reflectance difference 
A[py(X)]gg has been computed and is shown in figure 9. This figure indicates 
that the two background waters give differences in magnitude of A[p^(X)]gg, 
but its spectral variation is nearly independent of whether mixture 4 or Back 
River water is the background water. For this particular pollutant (Southwest 
Philadelphia Waste Treatment Facility sewage sludge) , the spectral variation 
of AtPu(X)]gg appears to be less sensitive to background-water variations than 
either fPu(X) 1 ss/^Pu^^^ ^ bw ^ ss/^Pu^^) ^ bw* 

To further investigate characteristics of the A[py(X)]gg parameter, com- 
putations were made to determine how spectral characteristics of A[p^(X)]gg 
might change with sewage-sludge concentration. These results are shown in fig- 
ure 10 for three types of background water based on data from figure 5 and ref- 
erence 9. The parameter A[pjj(X)]gg is shown for filtered, deionized tap 
water because laboratory reflectance measurements are not considered valid for 
clear background water, as discussed in appendix A. In clear water, laboratory 
spectral measurements are not usable unless a ^ 2 m“^ (concentration of 
suspended solids usually greater than or equal to 4 ppm) . In this case, 
A[Pu(X)lss used to illustrate changes which occur as sewage-sludge 
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concentration increases frcm 4 to 69 ppm, whereas A[py(X)]gg is used to study 
changes as sludge concentration increases from 0 to 69 ppm. Figure 10 shows a 
nearly uniform spectrum for both A[py(X)]gs and A[py(X)]gg for sludge 
concentrations between 4 and 69 ppm for all three types of background water. 

The small variations evident in some curves are within the standard deviation 
of laboratory reflectance values (o = ±0.16 percent). 

Knowledge of atmospheric effects is required to relate surface-reflectance 
to apparent-radiance parameters. Equation (3) shows that values of atmospheric 
transmission T^CX) and diffuse-reflector radiance Lgo(X) are required to 
convert the difference in inherent reflectance at the water surface A[pu(X)lgg 
to the difference in radiance at altitude A[L 2 (X)]gg. Typical values for 
TgCX) and LgQ(X) are shown in figure 11. Figure 11(a) is based on Lgo(X) 
measurements made on the beach at Cape Henelopen, Delaware, on July 24, 1978, at 
9:40 a.m. EDT. Original measurements were made with a spectrcmeter at a spec- 
tral resolution equal to 16 nm. The values shown in figure 11(a) are integrated 
average values which simulate each band of the M^s instrument. Values for 
Tg(X) were computed every 6 nm by use of the analytical model described in ref- 
erence 12. Likewise, the values shown in figure 11(b) are integrated averages 
which simulate the spectral resolution of the M^S instrument. Figures 11(a) 
and 11(b) show that Lgo(X) and Tg(X) typically have different spectral vari- 
ations. This means that A[L 2 (X)]gg will have different spectral characteris- 
tics than A[py(X)]gg, as shown in figure 12. Figure 12(a) shows laboratory 
values of A[pjj(X)]gg taken directly from figure 10(b). Calculated values 
of A[L 2 (X)]gg are shown in figure 12(b). Values for Lgo(X) and Tg(X) from 
figure 11 were used in equation (3) to obtain A[L 2 (X)]gg. The standard devia- 
tion shown for A[L 2 (X)]gg is based on the a value for A[p^(X)]gg multi- 
plied by Tg(X) and LgQ(X) . In the true case, both ®nd Lgo(X) 

also contain measurement errors which propagate into A[L 2 (X)]gg and increase 
the actual standard deviation beyond that shown in figure 12(b). From fig- 
ure 12, it is clear that calculated values of A[L 2 (X)]gg contain significant 
uncertainty because of the inaccuracy of basic radiance measurements and the 
number of analytical operations which must be performed to obtain quantitative 
results . 

It is of interest to compare A[L 2 (X)]gg values calculated from laboratory 
tests with results from an actual flight experiment. On October 7, 1976, an 
experiment was conducted off the coast of Delaware in which environmental con- 
ditions were similar to those used in the calculations for figure 12. An M^S 
instrument was flown at an altitude of 3.1 km over a dump of pr imary- treated 
sewage sludge from the Southwest Philadelphia Waste Treatment Facility. Time 
of overpass was 11:56 a.m. EDT and atmospheric visibility was approximately 
15 km. Solar elevation angle was 44°. The flight data from that test contain 
scane noise which is typical of most multispectral scanner instruments. The 
average standard deviation for bands 1 through 9 is ±0.115 x 10“^ mW/cm^-sr-nm 
for individual pixels representing an area of 9 m by 9 m each. To reduce this 
noise, new radiance values were calculated for each pixel based on the average 
of 49 surrounding pixels (a 7-pixel by 7-pixel square) . The pixel averaging 
process reduced the average noise standard deviation from ±0.115 x 10~3 mW/cm^- 
sr-nm to ±0.058 x 10“^ mW/cm^-sr-nm. The noise level differed significantly 
between individual bands. 
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To obtain A[l2(^)1ss' "averaged" radiance of a pixel in the sewage-sludge 
plume was subtracted from the "average" radiance of the background water. The 

standard deviation of A[l 2 {A.) 1 ss estimated as the square root of 2 times 

the standard deviation of background-water noise for each band. The result of 
these calculations is shown in figure 13(a). It is not known what the sewage- 

sludge concentration was in the plume pixel used in these calculations. Fig- 

ure 13(a) indicates that the uncertainty in the flight data is of the same 
order of magnitude as that in the A[l 2 (^) 1 ss values calculated from labora- 
tory data (fig. 13(b)). Spectral variation of A[L 2 (X)]gs from flight data is 
similar in both magnitude and shape to the laboratory curves for the lower con- 
centrations of sewage sludge. Inaccuracies in both sets of data make it such 
that an accurate comparison with firm conclusions is not possible. 


CONCLUDING REMARKS 

In an effort to improve understanding of the effects of variations in back- 
ground water on reflectance spectra, laboratory tests have been conducted with 
various concentrations of sewage sludge diluted with several types of water. 
Spectral-reflectance characteristics have been obtained and applied to selected 
data-analysis algorithms which may have potential for pollutant identification. 
From these results, the following conclusions are made: 

1. Reflectance spectra for sewage-sludge mixtures in turbid waters are 
dependent upon the reflectance of the background water. 

2. The ratio of sewage-sludge reflectance to background-water reflectance 
showed some spectral variation for different turbid background waters. 

3. Spectral differences also occurred in the ratio of "delta" reflectance 
A[py(X)]gg to background-water reflectance as background water changed. For 
this reason, caution is suggested when the ratio of the difference in sewage- 
sludge apparent radiance and background-water apparent radiance to background- 
water apparent radiance from aircraft or satellite remote-sensing data is used. 

4. The difference in reflectance between sewage sludge and background 
waters did not show major spectral variations with different turbid background 
waters. Improved knowledge of atmospheric effects is required before labora- 
tory measurements of inherent reflectance can be accurately related to apparent 
radiance measured from aircraft or satellite remote sensors. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 12, 1979 
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APPENDIX A 


LABORATORY AND EQUIPMENT 

The cylindrical steel water tank has a diameter of 2.5 m and a depth of 
3 m. The botton is concave as illustrated in figure 1. The tank interior 
is coated with a black phenolic paint that absorbs 97 percent of incident 
radiation over the spectral range of the measurements (400 to 980 nm) . For 
these experiments, the tank was filled to within 0.3 m of the top with about 
n 600 liters of water. 

The circulation system was designed to maintain a vertical and horizontal 
homogeneous mixture in the tank and to maintain in suspension particles up to 
about 70 ym in diameter (specific gravity of 2.6) . This particle size corre- 
sponds to fine sand. In order to accomplish these design goals, water is 
pumped from the drain at the bottom of the tank into a system of pipes which 
returns the water to the tank through two vertical pipes on opposite sides of 
the tank. The pipes empty just above the concave bottom. Water entering the 
tank through these pipes washes over the concave bottom, meets at a location 
away from the drain, and wells upward. Tests using tracer techniques and 
transmission measurements have confirmed that this circulation system provides 
a nearly uniform homogeneous mixture throughout the tank. For pollutants such 
as sewage sludge with specific gravities less than 2.6, the present laboratory 
setup can suspend particles larger than 70 ym in diameter. 

The filtration and deionization system includes a commercial fiber 
swimming-pool filter, an activated carbon filter, and a charged resin 
deionizer. These units were placed in water lines parallel to the main cir- 
culation system water lines and can be used separately or in any combination 
by using valves. The two filters remove particulates and dissolved organic 
materials from the water before it reaches the deionizer, where dissolved ionic 
substances are removed. After tap water is conditioned through this system, it 
contains less than 0.5 ppm of suspended solids and less than 2 ppm of dissolved 
substances . 

The light source is a solar-radiation simulator designed to approximate 
the spectral content of the Sun's rays. The radiation spectrum is produced by 
a 2.5-kW xenon short-arc lamp and transferred to the target plane through an 
optical arrangement inside the simulator and a collimating lens accessory. 

With the collimating lens accessory, the projected beam is collimated to a 
0.1 5-m diameter at 0.3 m from the simulator and has a ±2.5° collimation angle. 
For these experiments, the simulator was located approximately 6.0m from the 
water tank (fig. 1). At this distance from the simulator, the beam is about 
0.57 m in diameter. A mirror positioned 1.52 m above the water tank reflects 
the center of the beam to the water surface. The incidence angle with the 
water surface is 13° to avoid specular reflectance. The mirror is a first- 
surface mirror coated with aluminum and protected by an overcoat of silicon 
monoxide. It has a 0.3-m diameter and reflects an elliptical spot on the water 
surface which has a maximum diameter of 0.35 m. The simulator spectral input 
to the water surface is similar to, but not a precise duplicate of, sea-level 
standard solar-radiation curves often used in engineering calculations 
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(ref. 13). Figure 2(a) shows that the standard sea-level curves are rather 
variable, depending on the solar elevation angle. Figure 2(b) shows the simu- 
lator spectrum normalized to the solar spectrum at 600 nm and a solar elevation 
angle of 30°. These curves suggest that when laboratory measurements are made 
at a 32-nm spectral resolution, the input spectrum and possibly the output mea- 
surements are similar to those that would be expected in the field for a solar 
elevation angle on the order of 30°. The total intensity of the light hitting 
the water surface is approximately 8 percent of that in actual field conditions. 

The rapid- scanning spectrometer system consists of a spectrometer unit 
with a telephoto lens attachment and a plug-in unit with an oscilloscope and 
camera attachment. The spectrometer unit with telephoto lens attachment is 
mounted 2.43 m above the surface of the water (fig. 1). The spectrometer uses 
a Czerny-Turner grating monochromator without an exit slit. The spectral out- 
put of the monochrcmator is focused on the target of a vidicon tube, where the 
spectrum is stored as an electrical charge image. An electron beam periodi- 
cally scans the vidicon target to convert the charge image into an electronic 
signal. This signal is processed by the plug-in unit, which also functions as 
a controller between the spectrcmeter and the oscilloscope. The signal is dis- 
played on the oscilloscope and is photographically recorded with a camera. The 
spectrcmeter is designed to measure power per spectral bandwidth (spectral 
power) . The oscilloscope screen is used to show displacement of the instrument 
measurement. Oscilloscope displacement is proportional to spectral power as 
shown in the equation 

P(X) 

D = (AT) 

S 


The signal is internally processed in such a manner that the vertical-scale 
sensitivity factor S has a constant value over the wavelength range from 400 
to 980 nm. Values of S were obtained by the manufacturer using calibration 
procedures described in reference 14. (After receipt of the instrument, the 
manufacturer's calibration was checked in an approximate manner prior to the 
tests described herein.) The upwelled spectral radiance Lu(X) is defined as 

P(X) 

^u(^) = — (^2) 

where A is the area of the spectrometer entrance slit and 0 is the accep- 
tance solid angle of the spectrometer. Radiance values given herein are based 
on power received at the detector and are not corrected for losses through the 
telephoto lens. Tests with and without the lens indicate that such losses are 
much less than 5 percent for wavelengths between 400 and 980 nm. 

Conbining equations (Al) and (A2) results in 
DS 

Lu(^) = "TT (A3) 
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or 

D 

L„(X) = - (A4) 

K 

where 

Afl 

K = 

S 

Thus, upwelled spectral radiance is determined from oscilloscope displacement 
and the proportionality constant K, which is a function of the calibration 
factor S (which includes optical transmissivity) as well as acceptance 
angle and slit area A. 

In order to obtain spectral reflectance Pu(X), radiance measurements 
are made of a 99-percent diffuse horizontal reflector near the water surface 
ligo(X) . Values obtained are proportional to the spectrum being input to the 
water by the solar simulator. Reflectance is then computed by using the fol- 
lowing equation: 


Pu(>^) 


[Lu(X) ] 

bso ( 


(A5) 


By adjusting the slit area, spectral resolution of the spectrometer may be 
changed. For the tests described herein, all laboratory measurements were made 
with a spectral resolution of 32 nm. 

The instrument has been observed to experience daily variations in the 
calibration factor K, which affects absolute accuracy. According to instru- 
ment specifications, absolute accuracy of the measurements is believed to be 
±20 percent in the 400- to 600-nm range and ±12 percent in the 600- to 900-nm 
range. (A comparison of results from a number of laboratory tests at Langley 
Research Center tends to verify the manufacturer's specifications.) Included 
in the absolute error is a repeatability uncertainty of ±13 percent in the 400- 
to 600-nm range and ±3.5 percent in the 600- to 900-nm range. Discussions with 
the manufacturer indicate that these values are believed to be representative 
of 3a error bands. Because of these absolute errors, spectral radiances from 
tests conducted on different days usually differ somewhat in magnitude. How- 
ever, the overall shape of the relative spectrum over the wavelength range is 
consistent between tests conducted on different days. 

Results from laboratory tests are always limited, because the natural 
environment is not precisely duplicated. The effects of diffuse skylight, for 
instance, are not simulated in the laboratory used in this study. (The percent 
value of diffuse skylight to parallel sunlight varies from day to day in the 
natural environment.) Calculations with quasi -single- scatter ing and multiple- 
scattering optical models (ref. 15) have indicated less than 1 percent differ- 
ence in inherent reflectance between Sun-only and sky-only cases, however. For 

11 


II 



APPENDIX A 


this reason, it is believed that inherent-reflectance spectra measured in the 
laboratory are characteristic of those in the natural environment for the same 
water mixture. 

Measurements of upwelled radiance for clear waters in this laboratory are 
not considered valid. Data tabulated in reference 16 suggest that beam atten- 
uation coefficient ot should have a minimum value of 2.0 m“^ before bottom 
reflection can be ignored when water depth equals 3.0 m. Reference 17 examined 
the effects of solar spot size and concluded that a should have minimum val- 
ues between 2.0 and 4.0 m“^ to eliminate losses in signal when spot diameter is 
0.3 m. The only data that should be considered valid from this laboratory are 
those for turbid water mixtures with a greater than 2.0 m“^ . 
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RELATION OF ALTITUDE RADIANCE TO SURFACE REFLECTANCE 

From reference 18, reflectance at sea level and upwelled radiance at an 
altitude are related as follows: 

Lz(>^) = [Ta(X)]{[pu(X)] [Lso(X)] + Lrd(X) + Li-s(X)} + La(X) (B1 ) 

where 

L 2 (X) apparent upwelled radiance at altitude z, itW/cm^-sr-nm 

Ta(X) atmospheric transmission fraction 

Pu(X) inherent upwelled spectral reflectance above water surface, percent 

of input (relative to a 100-percent diffuse reflector) 

Ijso(X) upwelled radiance above water surface frcan a 100-percent diffuse 
reflector, mW/cm^-sr-nm 

Lrd(X) upwelled radiance from specular reflection of diffuse skylight, 
mW/cm^-sr-nm 

L^g(X) upwelled radiance from specular reflection of sunlight, mW/cm^-sr-nm 

La(X) upwelled radiance from light scattered into receiver by the 

atmosphere, nW/cm^-sr-nm 

Atmospheric transmission Tg(X) is a function of atmospheric constituents, 
altitude, and remote-sensor pointing angle (ref. 14); Lg(X) is a function 
not only of atmospheric constituents and solar elevation angle, but also of 
target and background reflectance (ref. 19); Lgo(X) is a function of solar 
elevation angle and atmospheric conditions; Li-d(X) is a function of atmo- 
sfiieric conditions and remote-sensor pointing angle; and Lj;s(X) is a func- 
tion of sea state and windspeed in addition to atmospheric and pointing-angle 
effects (ref. 18) . Thus upwelled radiance measured at altitude is a function 
not only of inherent reflectance of the water but also of altitude, remote- 
sensor pointing angle, atmosj^eric constituents, solar elevation angle, sea 
state, and windspeed. 

If the target pollutant and background water are located sufficiently 
close together and no differences in surface films exist, such that Tg(X), 
I*rd{^)f l^rs^^)' ba(X) are equal for both signals, an expression for the 

ratio of sewage-sludge radiance to background-water radiance can be written 

[Lg ( X) ] gg [Tq ( X ) ] { [ Q ( X) ] gg [ LgQ { X ) ] + Lj-(g(X) + Lj-g(X)T + Lg(X) 

= — (B2) 

[Lz(^)lbw fTa(X)]{[pu(X)]bwfbgo(X)] + Lj-d(X) + L^gCX)} + La(X) 
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Rearranging gives 
[Lz(X)1 

SS fPu( ^ ss/ 

(B3) 

B + C 

where 

B = [Lu(A)]b„ = fPuO^)lbwLso(^) (B4) 


La(>^) 

C = Li-d(^) + Lrs(^) + (B5) 

Ta(X) 


Thus, the ratio of sewage- sludge radiance to background-water apparent radiance 
is a function of altitude, remote-sensor pointing angle, atmospheric constitu- 
ents, solar elevation angle, sea state, and windspeed, as well as water reflec- 
tance ratio. 


Another algorithm which has seme potential for identification of marine 
pollutants is one in which the difference between the sewage-sludge radiance 
and background-water apparent radiance is expressed as a ratio to the apparent 
radiance of the background water. From equation (B1 ) this ratio can be 
expressed as 


A [L^ (X) ] gg 
fi'z ( ^ bw 


{A[Py(X) ] gg/ [ Pu ( X) ] } 


c 

1 + - 
B 


where 


A [L^ ( X) ] gg - [Lg (X) ] gg 
A[pu(X)]gg = tp|j(X)]gg 


[Lz (X) ] 
tPu(>^) Ibw 


(B6) 


(B7) 

(B8) 


This algorithm is similar to that of equation (B3) in that it is a function of 
altitude, remote-sensor pointing angle, atmospheric constituents, solar eleva- 
tion angle, sea state, and windspeed, as well as inherent-reflectance charac- 
teristics of the water. 

A third algorithm often considered for identification of pollutants is the 
simple difference between the sewage-sludge and background-water apparent radi- 
ances as defined in equation (B7) . Again assuming that the sewage sludge and 
background water are located sufficiently close together and no differences in 
surface films exist such that Tg(X), Lj-( 3 (X), Lj-g ( X) , and Lg(X) are equal 
for both signals, the algorithm for apparent-radiance difference can be written 
from equation (Bl) as 
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RELATION OF ALTITUDE RADIANCE TO SURFACE REFLECTANCE 

Frcxn reference 18, reflectance at sea level and upwelled radiance at an 
altitude are related as follows; 

Lz(X) = [Ta(X)]{[pu(X)] [Lgo(X)] + Lrd(X) + Lrg(X)} + La(X) (Bl) 

where 

Lz(X) apparent upwelled radiance at altitude z, iriW/cm^-sr-nm 

Ta(X) atmospheric transmission fraction 

Pa(X) inherent upwelled spectral reflectance above water surface, percent 

of input (relative to a 100-percent diffuse reflector) 

Lgo(X) upwelled radiance above water surface frcan a 100-percent diffuse 
reflector, mW/cm^-sr-nm 

Lrd(X) upwelled radiance from specular reflection of diffuse skylight, 
mW/cm^-sr-nm 

Lj-g(X) upwelled radiance from specular reflection of sunlight, mW/cm^-sr-nm 

Lg{X) upwelled radiance from light scattered into receiver by the 

atmosphere, iriW/cm^-sr-nm 

Atmospheric transmission Ta(X) is a function of atmospheric constituents, 
altitude, and remote-sensor pointing angle (ref. 14); Lg(X) is a function 
not only of atmospheric constituents and solar elevation angle, but also of 
target and background reflectance (ref. 19); Lgo(X) is a function of solar 
elevation angle and atmospheric conditions; is a function of atmo- 

spheric conditions and remote-sensor pointing angle; and Li-g(X) is a func- 
tion of sea state and windspeed in addition to atmospheric and pointing-angle 
effects (ref. 18). Thus upwelled radiance measured at altitude is a function 
not only of inherent reflectance of the water but also of altitude, remote- 
sensor pointing angle, atmospheric constituents, solar elevation angle, sea 
state, and windspeed. 

If the target pollutant and background water are located sufficiently 
close together and no differences in surface films exist, such that Tg(X), 
Lrd(^)» l'rs(^)' equal for both signals, an expression for the 

ratio of sewage-sludge radiance to background-water radiance can be written 

tLz(X)]gs [Ta(X)]{[Pu(X)]gg[Lgo(X)] + Lrd(X) + Lrg(X)> + Lg(X) 

= (B2) 

[Lz(X)]bw [Ta(X)]{[pu(X)]b„[Lgo(X)] + Lrd(X) + Lrg(X)} + La(X) 
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Rearranging gives 

[Ljj ( ^) ] ss ^ ^ ss/ ^ bw^® ^ 

B + C 


(B3) 


where 


B - [Ljj ( X) ] = [ Pq ( X) ] bw^so ^ 


(B4) 


La(^) 

C = Lrd(X) + Lj-s(X) + — (B5) 

Ta(X) 


Thus, the ratio of sewage-sludge radiance to background-water apparent radiance 
is a function of altitude, remote- sensor pointing angle, atmospheric constitu- 
ents, solar elevation angle, sea state, and windspeed, as well as water reflec- 
tance ratio. 


Another algorithm which has seme potential for identification of marine 
pollutants is one in which the difference between the sewage-sludge radiance 
and background- water apparent radiance is expressed as a ratio to the apparent 
radiance of the background water. From equation (B1 ) this ratio can be 
expressed as 


A [L 2 (X) ] 23 
^1*2 ( X) 


{ A [ Pjj ( X) ] ss/ ^ Pu ( ^ bw ^ 


C 

1 + - 
B 


where 


A [L 2 ( X) ] 22 - [L2(X)122 
A[Pu(X)] 2 s = fPu^^^^ss 


( X) ] bw 
fPu( ^bw 


(B6) 


(B7) 

(B 8 ) 


This algorithm is similar to that of equation (B3) in that it is a function of 
altitude, remote-sensor pointing angle, atmospheric constituents, solar eleva- 
tion angle, sea state, and windspeed, as well as inherent-reflectance charac- 
teristics of the water. 

A third algorithm often considered for identification of pollutants is the 
simple difference between the sewage-sludge and background-water apparent radi- 
ances as defined in equation (B7) . Again assuming that the sewage sludge and 
background water are located sufficiently close together and no differences in 
surface films exist such that T^CX), Lp^(X), Lj .3 ( X) , and ^^(X) are equal 
for both signals, the algorithm for apparent-radiance difference can be written 
from equation (B1 ) as 
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A[Lz(X)]gs = Ta(X)Lso(X) A[pu(X)]gs (B9) 

This algorithm appears to have some advantage over those of equations (B3) 
and (B6) in that variations due to specular reflection of skylight and sunlight 
as well as atmospheric backscattering have been eliminated. This means that 
radiance difference is a function of inherent-reflectance difference, altitude, 
remote-sensor pointing angle, atmospheric constituents, and solar elevation 
angle. 
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TABIiE I.- COMPONENTS OF ARTIFICIAL SEAWATER 


Sequence of 


Amount, grams, added 

addition 

Compound 

per liter of filtered. 



deionized tap water 

1 

NaCfi- 

12.265 

2 

Na 2 S 04 

2.045 

3 

kcA 

.350 

4 

NaHC 03 

.100 

5 

KBr 

.050 

6 

H 3 BO 3 

.015 

7 

MgCS-- 6 H 20 

5.550 

8 

SrCil 2 ’ 6 H 20 

.255 

9 

CaCfi-2 

.580 

10 

NaF 

.0015 


TABLE II.- NUTRIENTS ADDED TO ARTIFICIAL SEAWATER 


Compound 

Amount, grams, per liter 


of artificial seawater 

KNO 3 

0.00600 

Na 2 HP 04 

.01000 

FeS04*7H20 

.00498 

ZnS04*7H20 

.00882 

MnCS-2'4H20 

.00144 

M 0 O 3 

.00071 

Co(N 03 ) 2 * 6 H 20 

.00049 

EDTA 

.05000 

KOH 

.03100 

H 3 BO 3 

.01140 

Fe(NH 4 ) 2 ( 304 ) 2 • 6 H 2 O 

.03510 

NaEDTA 

.03300 

Biotin 

.00010 

B 12 

.00010 

Thiamine HCS- 

.02000 

Na2Si03’9H20 

.00466 
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TABLE III.- RESULTS OF LABORATORY ANALYSIS OF TEST-WATER MIXTURES 


Mixture® 

Description 

m ' 

Suspended solids 

Chlorophyll a 
concentration, 
Ug/il 

Algae 
count , 
cell s/m S, 

Organic, 

ppm 

Inorganic, 

ppm 

Total, 

ppm 

1 

Artificial seawater plus algae 

1.5 

3.6 

3.4 

7.0 

7 

• 

00 

X 

o 

2 

Mixture 1 plus algae 

2.2 

4.7 

2.1 

6.8 

15 

3.4 

3 

Mixture 2 plus algae 

2.9 

4.6 

5.0 

9.6 

20 

4.5 

4 

Mixture 3 plus overnight growth 

7.3 

6.2 

5.4 

11.6 

31 

6.0 

5 

Mixture 4 plus sewage sludge 

9.2 

8.9 

10.4 

19.3 

26 

6.5 

6 

Mixture 5 plus sewage sludge 

21 .2 

20.7 

20.0 

48.7 

27 

6.2 


^Mixture 0 is considered to be artificial seawater plus nutrients. 









Figure 1 Sketch of laboratory setup. 
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Solar elevation, deg 



(a) Standard sea-level solar irradiance spectra (ref. 13). 


Solar simulator (32-nm spectral resolution), past tests 

Solar simulator (32-nm spectral resolution), present tests 

Solar elevation. 30° 



Wavelength. X, nm 


(b) Solar-simulator and standard sea-level spectra. 

Figure 2.- Standard sea-level solar irradiance spectra and comparison with 

solar-simulator data. 
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Figure 3.- Laboratory test sequence. 
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(a) Spectral radiance. 
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(b) Spectral reflectance. 

Figure 4.- Laboratory spectral radiance and reflectance for artificial 
seawater-algae mixtures. Spectral resolution, 32 nm. 
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Sewage 

Chlorophyll, sludge, Total, 
Mixture pg/I ppm ppm 

4 31 0 11.6 



X, nm 


(a) Spectral radiance. 

Figure 5.- Laboratory spectral radiance and reflectance for background water 
and sewage-sludge mixtures. Spectral resolution, 32 nm. 
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Suspended solids 


P,(A). 

percent 


Chlorophyll^, 
Mixture |jg/ i 

-4 31 


Sewage 

sludge. 

ppm 

0 


Total. 

ppm 

11.6 



X, nm 

(b) Spectral reflectance. 
Figure 5.- Concluded. 
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Fig. 5 data (mixture 6) 



(a) Sewage- sludge mixtures of about 34 ppm. 


Fig. 5 data (mixture 4) 



(b) Background water. 

Figure 6.- Spectral reflectance of mixture 4, mixture 6, and Back River water. 

Spectral resolution, 32 nm. 
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(a) Laboratory-measured spectra. 
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(b) Simulated spectra for M^S instrument. 



1000 


7.- Ratio of sewage-sludge reflectance to background-water reflectance 
for sewage-sludge concentration of about 34 ppn. 
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(a) Laboratory-measured spectra. Spectral resolution, 32 nm. 



Figure 8.- Ratio of difference between sewage-sludge reflectance and 
background- water reflectance to background-water reflectance for 
sewage-sludge concentration of about 34 ppm. 
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(a) Laboratory-measured spectra. 
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(b) Simulated spectra for M^S instrument. 

Figure 9.- Difference between sewage-sludge reflectance and background-water 
reflectance for sewage-sludge concentration of about 34 ppm. 
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(a) Mixture 4. 
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(b) Back River water. 
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(c) Filtered, deionized tap water. 

Figure 10.- Spectra of reflectance differences for several sewage-sludge 
concentrations in different background waters. Spectral resolution, 
32 nm. 
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(a) Radiance from diffuse reflector measured on July 24, 1978, 



(b) Calculated atmospheric transmission frcan analytic model of reference 12. 

Figure 11.- Typical diffuse-reflector radiance and atmospheric transmission 

spectra for M^S instrument bands. 
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(a) Reflectance differences for Back River water laboratory tests. 
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(b) Calculated radiance differences. 

Figure 12.- Ccsmparison of reflectance differences and radiance differences 
for M^S instrument and assumed atmospheric conditions. 
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(a) Aircraft M^S instrument data. 
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(b) Calculated values from laboratory results. 


Figure 13.- Characteristics of field and extrapolated laboratory results. 
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